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Abstract. In this paper, a dc Josephson junction between two singlet superconductors (d-wave and 
s-wave) with arbitrary reflection coefficient has been investigated theoretically following the famous 
paper [ Y. Tanaka and S. Kashiwaya 1996 Phys. Rev. B 53, R11957]. For the case of High T c 
superconductors, the c-axes are parallel to an interface with finite transparency and their ab-planes 
have a mis-orientation. The effect of transparency and mis-orientation on the currents is studied both 
analytically and numerically. It is observed that, the current phase relations are totally different from 
the case of ideal transparent Josephson junctions between d-wave superconductors and two s-wave 
superconductors. This apparatus can be used to demonstrate d-wave order parameter in High T c 
superconductors. 



PACS numbers: 74.50.+r,85. 25. Cp,74.20.Rp,85.25.Dq, 74.78. Bz 

The quantum transport through the junctions between high T c unconventional superconductors has 
attracted much attention in recent years, particularly in view of various recently discovered copper 
oxide compounds. It is well-known that the isotope effect for these superconductors is rejected or weak. 
The isotope effect is the most important reason for the phonon participation in the superconductivity 
mechanism. Thus we should search for another pairing mechanism for high-T c superconductors. In 
spite of enormous efforts made to understand the physical mechanisms of pairing in various high- 
Tc superconductors, the situation still remains unclear. A central role in understanding of this type 
of superconductivity relates to the symmetry of the order parameter. Soon after the first Bednords 
and Muller [T] discovery, the symmetry of the d a; 2_ !/ 2-form was suggested for high-T c superconducting 
materials. It was considered that the formation of midgap states (MGS) at surfaces and interfaces 
of d-wave superconductors affects the current transport properties of junctions involving d-wave 
superconductors [2]. It is found that in spite of s-wave superconductors for <i-wave superconductors 
we can observe midgap states. This phenomenon is due to the changing the sign of order parameter in 
d-wave superconductors. On the other hand, measurements of direct Josephson current yield valuable 
information on the symmetry of the order parameter which is essential for understanding the mechanisms 
of superconductivity in these complex materials. Phase interference experiments definitely suggest the 
presence of d-wave symmetry of the order parameter in high-T c superconductors [H [4] . The properties 
of spontaneous currents generated at surfaces and interfaces of d-wave superconductors using the self- 
consistent quasiclassical Eilenberger equations are investigated [5]. The effect of transparency and 
mis-orientation on the current as the main purpose of this paper. In this paper, we have solved the 
nonlocal Eilenberger equations [6 supplemented with Zaitsev boundary conditions [7] and obtained 
the corresponding Green functions analytically. These boundary conditions have been rewritten in 
[El [9]. Then, using the Green functions, we compute the spontaneous and Josephson currents. The 
system which is investigated in this paper, consists of two d-wave superconductors connecting through 
an interface. The interface has the finite transparency for quasiparticles and two superconductors have 
a mis-orientation with each other. Quasiclassical theory of superconductivity is based on the nonlocal 
Eilenberger's equations for the quasiclassical matrix propagator. In the case of a clean (ballistic) singlet 
anisotropically paired superconductor, the equations for the ^-integrated Green functions reduce to the 
following 2x2 matrix form [6] 

d ~ - 
vf • ^-Guj{vf, r) + [ljt 3 + A(v F , r), G w (v F , r) = (1) 
or 
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Figure 1. Josephson current through the interface versus phase for D = 0.99, on = 0, T/T c = 0.1, 
jo = 4neN(0)vpT c . 



where 



A 
At 



G u (v F ,r) 



9w 



fuJ 



(2) 



A is the superconducting order parameter, rgis the Pauli matrix, and G w (vf,i") is the matrix Green 
function which depends on the electron velocity on the Fermi surface v>, the coordinate r and the 
Matsubara frequency to — (2n+l)irT ,with n and T being an integer number and temperature respectively. 

We also need to satisfy the normalization condition, = \Jl — fufui with being time-reversal 
counterpart of f u . In general, A depends on the direction of v F and is determined by the self-consistency 
equation 



A(v F ,r) = 2nN(0)T V (V(v F , V F )Uv' F , r)> , 



LJ>0 



(3) 



where V(vf,v' f ) is the interaction potential. Solution of the matrix equation ([T]) together with self- 
consistency equation (|3|) and normalization condition determines the current density j(r) in the system 

j(r) - -47rzeiV(0)T^ (v F g u (w F , r)) vp . (4) 



2tt 

In two dimensions, N(0) — m/27r is the 2D density of states and (...) = J (d6 /2tt)... is the averaging 

o 

over directions of 2D vector v F . A symmetric configuration is considered and the equations are changed 
to a two-dimensional problem. Thus the spatial distributions of A(r) and j(r) depend only on the 
coordinates in the a6-plane, and Eilenberger equations |T]) reduce to a simpler problem. In the present 
paper, it is assumed that transparency of the interface to be finite 0<J5<1. Also, the equations ([1]) for 
Green function G w (v^,r) have to be supplemented by the suitable boundary conditions [7j|8l[9]. These 
boundary conditions couple the transmitted and reflected trajectories [5]. So, because of finite reflection 
at the interface, it is impossible to impose the continuity of the Green functions along a trajectory that 
crosses over the interface, as done for the ideal transparent system |17) . Instead, the Zaitsev boundary 
conditions [7] should be exerted as follows 
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Figure 2. Parallel, spontaneous current j y (0) versus phase <j> for T/T c = 0.1, a; = 0, D = 0.99 and 
different a r s. 



where 

s r = GJ(v F ,x=0)+GJ(v' F ,x=0) (7) 

and 

d r = GJ(v F ,x=0)-GJ(v' F ,x=0) (8) 
with being the reflection of v^with respect to the boundary and D is transparency coefficient of the 
interface that can be momentum dependent generally. Similar relations also hold for s and d 1 . In 
our calculations, Ai yr (v F ) — Ao(T)cos(2(# — a l,r)) are used for the left and right side of the interface 
respectively with d x 2_ y 2 order parameters symmetry and a;. r being orientations with respect to the 
interface. It is remarkable that a step function for the spatial dependence for the phase of order parameters 
are considered as A^ r (x) — A; jr exp ±i(f>/2 respectively and for simplicity the self-consistency of problem 
is ignored. For all quasiparticle trajectories the Green functions satisfy the boundary conditions both 
in the right and left bulks and at the interface. We have solved this problem and obtained the Green 
functions of this system analytically. Using our solution of Eilenberger equations and symmetries of 
singlet order parameters, we obtain the following expression for the current density ([4]) exactly at the 
interface 

■ rm . Zsr/ v77A ; A r sin0 \ 

J ° T c ^ \^ + A^ r cos ct>+C-^nn r /^ W 



where, £\ r = yu> 2 + Af r , with u> being the Matsubara frequency, v = v F /v F is the unit vector, 
r) = sign(v x ), <fi is the phase difference between two bulks and jo = 47reiV (Q)v F T c . 

For the junction between two spin-singlet s-wave superconductors, with the isotropic order parameters 
in the momentum space A(v F ) = A, at the junction we reproduce the current density as follows: 




2 



K0) = ^r^_==tanh-Wl- jD ( S in5) . (10) 



j A Dsin 

^-^(sinf) 

Where, D is assumed to be independent of momentum direction and isotropic in the momentum space. 
This formulate coincides with the results of [HI [T31 HH Q~5] and for ideal transparent point contact with 
results of [THl HZ] . At the interface and near the critical and zero temperatures we have 
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Figure 3. Josephson current through the interface versus phase for D = 0.5, a; = 0, T/T c = 0.1 and 
different a r s. 



and 

J(r~o) = *ff 0) . D " in " , (12) 

respectively. On the other hand, for the small values of transparency of interface between d-waves, the 
current expression is the following 

m = jo Y Y,\ — — A (13) 

that is the same as one of the relations in [10]. As it is clear from (|9|), for high values of transparency, 
currents are nonlinear functions of the transparency, D, because of (^-gp) term in ©. While for the low 
values of transparency currents are exactly linear functions of transparency (|13p . On the other hand, 
when nonlinearity is removed, the deviation from the sinusoidal form disappears. It is well-known that 
the current in this system for the low temperatures does not have sinusoidal form, while near the critical 
temperature it tends to the sinus [17]. Also, this system at the low transparency, has a current-phase 
dependence similar to the high temperature. For the Josephson junction between d-wave superconductors, 
in addition to the Josephson current another term of current is observed, called the spontaneous current. 
The spontaneous current that is parallel to the interface, is plotted as a function of mis-orientation and 
transparency in (Fig[2j) and (FigH]). The critical value of the spontaneous current is plotted in terms of 
transparency in (Fig[5j) . The Josephson and spontaneous currents are plotted as the functions of <fi for 
different mis-orientations at T/T c = 0.1, D = 0.99 (FigfTJ), (FigfS]) and for a low value of transparency 
D = 0.5 in (FigO and (FigH]). By decreasing the transparency, the currents decrease (Fig(5j). For 
high values of transparency the currents are not sinusoidal function of phase, but by decreasing the 
transparency they tend to the sinus as is observed in relations ([9]), IT3]) . The period of current-phase 
plot is interesting. We fixed the left orientation at ai = 0, by increasing the right orientation from to 
?, period of the current-phase graphs varies from 2?t to ir then, if the right orientation changes to ^, 
this period returns to 2n (FigfTJ) and (FigfSJ). The critical values of Josephson current in terms of mis- 
orientation are periodic function of left and right mis-orientations. They are periodic function with the 
period of w, because cos(2(a; + tt)) = cos(2a;). For high values of temperature, T c — T<^T C , the Josephson 
current has been obtained analytically and precisely as 

4tt(1 - X) 

1/(0) = Jo 315C(3) C #sin(0)[15cos(2a; - 2a r ) - cos(2a ; + 2a r )]i (14) 
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Figure 4. Parallel, spontaneous current j y (Q) versus phase rfr for T/T c = 0.1, a; = 0, D = 0.5 and 
different a r s. 



and for the spontaneous current we have 
j 167r(l - £•) 

j»(0) = 3i5 C ( 3 ) ^sin0sin(2a t + 2a r ) (15) 

where C( n ) is the Riemann zeta function. For the d-wave order parameters at high temperatures, 
T c — T<CT C , an analytical relation A(T) = J ^ticfz) ~ ( ^ ~ T~)' nas been used. It is observed that close to 
the T c , the currents are exactly sinusoidal function of phase, quadratic function of order parameter (linear 
dependence of (T c — T)) and the linear function of transparency, (Fig(SJ), (p3J and (fTB")) . The currents 
dependence on the orientations cti and a r are trigonometric functions and consequently periodic with the 
period of it generally. It is obtained that, in addition to the mis-orientation \ai — a r \, the orientation with 
respect to the interface is important, specially for the spontaneous current because of the terms including 
\ai + a r \ in the currents expressions. The spontaneous current for the high values of transparency is an 
abnormal function of the phase, but for the small values of transparency, it gets the sinus form (FigJ2J) 
and (Fig®. 

In conclusion, the effect of mis-orientation and transparency on the currents through the Josephson 
junction in the D-I-D system as well as in the S-I-S as an example, are investigated. The authors of 
[TT] . argued that, j(0 + tt) = —j((p), from my Green function is not generally correct, because of cos(</>) 
term in the denominator, but for the high values of the temperature and small values of transparency it 
will be correct. For the fixed left orientation, the current is a periodic function of the right orientation 
and the vise versa. In addition to the mis-orientation, \ai — a r \, the angle with the interface plays a role 
in the current, because j(cn,cti) ^ + (3, cti + (3), This case is returned to the scattering from the 
interface. The effect of transparency on the currents, as in equations (fT"3")) . (fH)) . (fT5|) and my numerical 
calculations (Figj5]) , for the high values of temperature and low values of transparency is linear. But for 
high values of transparency and the low temperatures the effect of transparency is nonlinear Thus 
nonlincarity of currents in terms of transparency j(-D) is coupled with the nonsinusoidal form of currents 
versus phase difference between order parameters j(0). These currents dependence on the transparency, 
phase difference and mis-orientation can be applied to show the d— wave pairing symmetry for the high 
T c superconductors. As the most important result it is obtained that j(ai, a r , 4> + 7r) = j(<%h a r + f , 4>)i 
that means changing the sign of order parameter on the Fermi surface changes the zero junction to the 
7r junction. 
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Figure 5. Critical Josephson and spontaneous currents versus transparency for oj; = 0, a r = \ and 
T/T c = 0.1,0.9. 
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